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ABSTRACT 

Recently, iDe Pontieu et all ()200 7lj) discovered a class of spicules that evolves more rapidly than 
previously known spicules, with rapid apparent motions of 50-150 kms~^, thickness of a few 100 km, 
and lifetimes of order 10-60 seconds. These so-called type II spicules have been difficult to study 
because of limited spatio-temporal and thermal resolution. Here we use the IBIS instrument to search 
for the high velocities in the disk counterpart of type II spicules. We have detected rapidly evolving 
events, with lifetimes that are less than a minute and often equal to the cadence of the instrument (19 
sees) . These events are characterized by a Doppler shift that only appears in the blue wing of the Ca 
II IR line. Furthermore the spatial extent, lifetime, and location near network, all suggest a link to 
type II spicules. However, the magnitude of the measured Doppler velocity is significantly lower than 
the apparent motions seen at the limb. We use Monte Carlo simulations to show that this discrepancy 
can be explained by a forward model in which the visibility on the disk of the high-velocity flows 
in these events is limited by a combination of linc-of-sight projection and reduced opacity in upward 
propelled plasma, especially in rcconncction driven jets that are powered by a roughly constant energy 
supply. 

Subject headings: Sun: chromosphere — Sun: atmospheric motions 



1. INTRODUCTION 

The solar chromosphere, when viewed at the limb, is 
dominated by a variety of jet-like features that have 
been difficult to understand, mostly because the spatio- 
temporal resolution of previous observations has been in- 
sufficient to resolve much of the dynamics, w hich has 
led t o a multitude of different interpretations (|Sterlind 
I2OOOI ). Recently observations at high resolution (100-200 
km) and cade nce (a few seconds) at the Swedish 1-m So- 
lar Telescope (|van Noort fc Rouppe vari der VoortJ[200fih 
and wit h the Solar Optical Te lescope (SOT) onboard 
Hinode ()De Pontieu et aO2007b[ ) have provided unprece- 
dented views of the complex mix of highly dynamic and 
finely structured features that dominate the chromo- 
sphere. A full understanding of chromospheric energetics 
and coupling to the transition region and corona requires 
a detailed study of what drives this plethora of features 
and how they are connected to one another. 

Some progress in this direction has been made by 
iDe Pontieu et al l ()2007bD who found that there are two 
fundamentally different types of spicules at the limb, 
with very different dynamic properties. The type I 
spicules evolves on timescales of order several minutes 
and seems to be driven by shock waves that form when 
oscillations and ffows leak into the chromosphere along 
magnetic flux concentrations. These spicules seem to 
correspond to dynamic fibrils (active regions) and a 
subset o f mottles (quiet Sun) when seen on the so- 
lar disk (jHansteen et al.l 120061 : iDe Pontieu et al.ll2007at 

Electronic address: oysteol@astro.uio.no 

^ Institute of Theoretical Astrophysics, University of Oslo, P.O. 
Box 1029 Bhndern, N-0315 Oslo, Norway 

^ Lockheed Martin Solar and Astrophysics Lab, 3251 Hanover 
St., Org. ADBS, Bldg. 252, Palo Alto, CA 94304, USA 

^ Center of Mathematics for Applications, University of Oslo, 
P.O. Box 1053 Blindern, N-0316 Oslo, Norway 

INAF-Osservatiorio Astrofisico di Arcetri, 1-50125 Firenze, 
Italy 



iRouppe van der Voort et al.l 120071 ). Hinode/SOT obser- 
vations revealed a second class of spicules ("type H") 
that occur on timescales of just a few tens of seconds 
during which they are seen to rise and then rapidly dis- 
appear, possibly because they are heated out of the band- 
pass. These newly discovered features, which dominate 
the li mb in coronal holes and ( to a lesser extent) the quiet 
Sun ()De Pontieu et al.ll2007bD . remain mysterious since 
the line-of-sight projection makes it difficult to determine 
whether the high apparent velocities (50-150 kms~^) are 
associated with real plasma flows of the same magnitude. 
There is a signiflcant literature on the properties and re- 
lationship of mottles and spicules in various wavelengths. 
However, the enormous advances in spatial and temporal 
resolution of the Hinode data, and the lack of high ve- 
locity observations in the classical observations limit the 
usefulness of detailed comparisons with previous obser- 
vational reports. 

In this letter we present a candidate for the disk coun- 
terpart of type II spicules. We observe short lived blue- 
shifted excursions in the line proflle of the Ca II 854.2 nm 
line that are not succeeded by an obvious re dshift. Sim- 
ilar events have been reported earher (e.g IWang et al.l 
ll998t[Chae et al.lll998f) . but on temporally and spatially 
more extended features. To explain the differences ob- 
served between the limb and the disk we present numeri- 
cal experiments using radiative transfer calculations and 
Monte Carlo simulations. 

2. OBSERVATIONS 

The observations were obtained using the Interferomet- 
ric BIdimensional Spectrometer (IBIS) mounted at the 
Dunn Solar Telescope (DST) of the National Solar Ob- 
servatory. For a mo r e tho rough description of the instru- 
ment see C avallin ii (2006'). The data used in this letter 
was obtained on 2 June 2004, during a period of good 
to excellent seeing conditions lasting for about 55 min- 
utes (UT 14:59-15:54). For a detailed description of the 
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Fig. 1. — A snapshot of a line center image showing an overview 
of the field of view. Crosses indicate the positions of the observed 
rapid blueshifted excursions, with the five examples seen in Fig[2] 
marked as diamonds with the corresponding number. 

data set and the reduction routines we ref er the reader 
to ot her work using the same data set (e.g. iCauzzi et al.l 
120081 ). The data was obtained with a scanning sequence 
of 27 steps in the Ca II 854.2 nm hne, with a step size 
of 8-16 pm and a spectral transmission FWHM of 4.4 
pm. This scan covers a spectral region of -0.12 nm to 
about 0.16 nm from line center. The signal-to-noise ra- 
tio is about 130 in the linewing. The diffraction limit 
of the telescope at 854.2 nm is ^ = 0.23" and the pixel 
size is 0.082" but with 2x2 binning the effective size is 
0.165". The scanning of the Ca line takes 7 seconds and 
was repeated every 19 seconds pointing at an enhanced 
network very close to the disk center. The field of view 
(FOV) covers four rosette structures (see FiglT]). 

3. DATA ANALYSIS 

Visual inspection of the data set shows that the general 
dynamic behavior of the spectral line varies dramatically 
across the different regions in the FOV (e.g. network , 
internetwork (IN), and canopy, see iVecchio et al][20?)l . 
The IN is cluttered with brightenings caused by th e 3 min 
shocks as explained by ICarlsson fc Steiril (|1997[ ) , while 
the center of the rosettes are covered with mottles (with 
a sp ectral signature very s imilar to that of dynamic fib- 
rils, |Langangen^etIal] [200^. At the edges of the rosettes 
we observe a marked difference in the dynamic behavior 
of the line profile. The line is wider and there are no 
dominating shocks, as compared to the network and the 
IN. It is in this region, at the edge of the rosettes that 
we find the blue shifted excursions (see Fig[T]and Figl2|). 

The events have been identified in a semi-automatic 
fashion. First we perform a search for changes in the in- 
tensity on the blue side of the line at around 20 kms^^. 
We then exclude all events that show a succeeding red- 
shift, since these kinds of events are identical to the tem- 
poral evolution of the velocity in a dynamic fibril or an 



IN shock. Finally we exclude all events that occur during 
times of worse seeing quality. Such an automatic search 
results in a manageable number of events. The events 
found by the search routine were manually inspected and 
only clearly blueshifted and short lived events were cho- 
sen for further analysis. In this fashion we identified 
87 blue shifted events distributed over the whole time 
series, from now on denoted as rapid blueshifted excur- 
sions (RBE). We have done a visual search for separate 
downflows associated with the RBEs and have found no 
clear candidates. 

To measure the spatial extent of the RBEs we calcu- 
late the spatial correlation of the dynamic behavior of 
the line profile at the semi-automatically flagged loca- 
tion with that of surrounding pixels. We choose the most 
pronounced timestep in the RBE as anchor, subtract the 
mean line profile from the spectra, and calculate the cor- 
relation coefficient of the temporal evolution of the blue 
side of the RBE between neighboring pixels. The area 
covered by pixels with higher than 80 % correlation to 
the anchor pixel is used as a measure of the spatial ex- 
tent of the structure. The RBEs can be elongated (78%), 
round (8%), or have an irregular shape (14%), see Figl21 
The length of the elongated RBEs is usually 0.5-1.5 Mm 
with an average of 1.2 Mm and the width is usually 0.3- 
0.6 Mm with an average of 0.5 Mm. The small round 
RBEs are often not more than 0.5 Mm in diameter. It 
is likely that some of the elongated shape of the RBEs is 
caused by the smearing due to the long scanning times 
(7 sees) of the data. 

From the measured RBEs we derive a mean lifetime 
of 45 ± 13 seconds. It is difficult to accurately measure 
plasma velocities in the RBEs, since many of the spec- 
tral profiles seem to consist of more than one compo- 
nent (but without any separation between components) . 
However, a rough estimate of the velocity based on visual 
inspection and comparison to radiative transfer calcula- 
tions suggests velocities of order 15-20 kms^^ (see 21for 
more discussion of this topic) . 

4. NUMERICAL EXPERIMENTS 

The similarities in properties such as uniquely upward 
motion, rapid disappearance, lifetime, and location sug- 
gest a link between type II spicules and RBEs. However, 
there is an apparent lack of very high velocity RBEs, es- 
pecially compared to the 50-150 km s~^ velocities seen 
in type II spicules. To investigate these differences in the 
observed properties of type II spicules and RBEs we per- 
form several numerical experiments. First we construct 
a toy model for the RBEs: a 100 km wide component 
with a LOS velocity of 15 and 20 km/s is introduced in a 
FAL C model (Fon tenla et al.lll993f l. at a variety of dif- 
ferent heights. We then use MULTI (|Carlssonlll98i) to 
calculate the impact of these higher velocity components 
on the Ca II 854.2 nm line profile (as seen on the disk). 
While this may be unrealistic, it can give some insight 
into the effect of atmospheric motions on the line profile. 
A full understanding of the radiative transfer in these 
events will re quire a more sophisticated approach e.g. a 
cloud model ([Heinzel fc Schmiedeilll99^ . or more likely 
non-LTE radiative transfer from a 3D MHD model. This 
numerical experiment indicates that events that occur at 
the top of the chromosphere (where lower densities oc- 
cur) are increasingly more difficult to observe. 
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Fig. 2. — Five examples of RBEs are shown. The top row shows the temporal evolution of the RBE (arrow). The detailed spectra 
(central row) of the RBEs (solid line) is shown together with the mean profile (dotted line). The line center image is shown in the bottom 
row. The contour (white line) shows the spatial extent of the RBEs. Note the alignment of the RBEs with the field lines, see also Fig[T] 



As a next step we investigate whether line-of-sight pro- 
jection is behind the lower velocities (15-20 km s""'^) ob- 
served in RBEs on the disk compared to the high veloc- 
ities of type II spicules at the limb. Perhaps the strong 
inclination from the vertical of the magnetic field lines 
at the edge of rosettes causes a reduction in observed 
velocities along the line-of-sight? For this purpose we 
perform a Monte Carlo simulation in which 50 type II 
spicules are allowed to occur (randomly) during a time 
period of 55 min (at 19 s cadence). Each spicule has a 
lifetime randomly chosen from a Gaussian distribution 
around 40 s with a standard deviation of 20 s. During 
their lifetime, the spicules are assigned an upward ve- 
locity chosen from a Gaussian distribution around 70 ± 
15 kms~^. The spicules have random orientation with 
a uniform distribution of angles between the spicule axis 
and local vertical between 0° and 90°. We then synthe- 
size the impact of these spicules on the Ca II 854.2 nm 
line by assuming each type II spicule is associated with 
a Gaussian absorption component with a depth of 10% 
of the continuum value, with a half- width of 10 kms~^. 
The result of such a simulation (Figl3] A) shows a large 
number of high velocity events with Dopplcr shifts larger 
than 30 km/s, which we do not see in the IBIS data. We 
should note that assigning spicules orientations that are 
more inclined from the vertical (e.g., as in the inclined 
regions at the edge of rosettes) does lead to a reduction 
of high-velocity events in the simulated IBIS Doppler- 
grams. However, Hinode/SOT data shows that many 



type II spicules are more vertical, so a skewed angle dis- 
tribution cannot be the reason for the mismatch in ve- 
locities between type II spicules and RBEs. 

Next we investigate whether the decreasing density 
with height in spicules could lead to a reduced visibility 
of nearly vertical spicules in which the density rapidly 
decreases with height. A simple way of incorporating 
this effect is to decrease the absorption caused by each 
spicule as a funct ion of height, using the o bserved inten- 
sity scale height (|De Pontieu et al.ll2007"bl ). To test the 
maximal effect we use the lower observed scale height of 
1500 km. The absorption component is then allowed to 
decrease in strength according to the fraction of the scan- 
ning time (7 sees) and the cadence (19 sees) the spicule 
occupies at the different heights. Highly inclined events 
(e.g., at the edge of the rosettes) do not reach high up, 
so do not suffer much attenuation of the default absorp- 
tion. An example of such a simulation (FigOB) shows 
that the visibility of events with high Doppler shifts is 
reduced significantly. However, it appears there are still 
a larger number of such events than we observe in the 
IBIS data. 

A resolution to this problem comes when we assume 
that the spicules are caused by a local release of en- 
ergy, such as in a reconnection event. If the reconnection 
events occur over a range of heights and each event is 
assumed to release a roughly constant amount of energy, 
then the dropoff in density with height will cause events 
that occur in the upper chromosphere to naturally lead 
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Fig. 3. — Results from the Monte Carlo simulations; a simulation with infinite scale height (panel A), a simulation with a scale height of 
1500 km (panel B), and a simulation with a scale height of 1500 km and random height of energy input panel C) are shown. The vertical 
dashed line is the limit of the IBIS observations. 



to spicules with lower density and higher velocity than 
events driven by reconnection lower down. We incorpo- 
rate this idea into the Monte Carlo simulations by allow- 
ing each spicule to be driven by reconnection at a ran- 
dom height chosen from a uniform distribution between 
1000-2000 km (in the FAL C model). The velocity of the 
spicule caused by the reconnection is then assumed to be 
given hy V = y^2Uk/ p, where Uk is the kinetic energy 
density in the event, and p is the density in the FAL C 
model at the height of the event. The energy density of 
the event is now a free parameter, but there are a whole 
range of values for this parameter (10-70 ergscm"'^) that 
produce simulated Dopplergrams that are very similar to 
the RBEs in the IBIS data. One example of such a sim- 
ulation with Uk = 25 crgscm"'^ is shown Figl3]C. While 
this simulation produces a whole range of high velocity 
events, these are not visible against a bright background, 
and only visible at the limb (as type II spicules) where 
there is no background radiation. This can explain the 
lack of high velocity RBE's on the disk, and the predom- 
inance of high velocity events at the limb (where there is 
a bias towards mostly vertical high velocity events that 
reach greater heights). 

5. CONCLUSIONS 

We have presented a candidate for the disk counter- 
part of type II spicules: rapid blue-shifted events (RBEs). 
The RBEs are found near the edges of the network, where 
movies show that straw-like features are present as well. 
The lifetimes of the RBEs are about 45 seconds, which is 
very close to those of type II spicules (40 s) . The length 



of the elongated RBEs is about 1.2 Mm and the width is 
0.5 Mm. The length of the RBEs is shorter than the type 
II spicules, which might be caused by the intrinsic differ- 
ences in visibility between the two kind of observations. 
Since RBEs are identified by a Dopplershift they are as- 
sociated with mass motion. The magnitude of the mass 
motion in RBE's is lower than the apparent motion ob- 
served in the Type II spicules. Monte Carlo simulations 
suggest that this discrepancy is to be expected if there 
is an inverse relationship between density and velocity of 
the events; high density, low velocity events show enough 
absorption to be visible in IBIS data taken on the disk, 
whereas low density, high velocity events are more clearly 
visible at the limb (since they reach greater heights), but 
have little visibility on the disk because they have little 
opacity. Such an inverse relationship would be expected 
if type II spicules are caused by reconnection events that 
occur over a range of heights in the chromosphere, and 
in which the amount of energy is roughly constant from 
event to event. Such a scenario is not unrealistic, and is 
supported by the presence of RBEs at the edge of rosettes 
where the relentless magnetoconvection pushes the ubiq- 
uitous IN flux into the strong network flux concentra- 
tions. Compounding factors that help explain the mis- 
match in observed velocities include the effects of vari- 
able LOS velocity due to the magnetic field topology and 
the fall in intensity due to the decrease in density as the 
spicules reaches greater heights. We show that a combi- 
nation of these factors can explain the visibility of RBEs 
on the disk and type II spicules at the limb. 
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